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Abstract. The model of a geometrically thin gaseous disk in the external gravitational potential is 
considered. The dinamics of small nonaxisymmetric perturbations in the plane of the accretion disk 
with dissipative effects is investigated. It is showed, that conditions of development and parameters of 
unstable oscillation modes in the opticaly thick accretion disk are strongly depended on the models of 
viscosity and opacity. 

The possibility of the development of varions types of instabilities in a thin gaseous 
disks is very attractive one for understanding different aspects of the accretion-disk (AD) 
phenomenon. In order to explain the required large values of dissipative coefficients, 
the concept of the turbulent viscosity has been used, which may be caused by the 
developed turbulence of a gaseous medium arising from the loss of stability. On the 
other hand, the nonlinear evolution of unstable oscillation modes may be responsible 
for many nonstationary phenomena in accreting systems. 

There are four unstable oscillation modes in the framework of the standard a-model 
of an accretion disk [1], Two of then are acoustic [6-11], one is viscous and the other one 
is thermal [2-11]. A distinctive feature of the dynamic viscosity rj = an in the model 
of the accretion disks is it’s dependence on the surfase density a and on the disk half- 
thicness h [n is the kinematic viscosity). The perturbation of the dynamic viscosity rj is 
responsible for the formation of all unstable four modes. However, in all the above-cited 
papers on the dynamic of linear perturbation is suggested that viscosity simultaneously 
changes with changing of accretion-disk parameters. 

In our present work, we extend our research to time delay of the viscosity i nfl uence 
on the thermal, viscous and acoustic instabilities. In the construction of a different 
viscous models AD is suggested that viscosity caused by the developed turbulence p ~ 
autit, where ut and £t is the characteristic velocity of the large-scale turbulent pulsations 
and it’s characteristic size, respectively [12], The fundamental energy is contained in 
the large-scale pulsations, however dissipation of the energy is the case in the small- 
scale pulsations. By this means, as the local conditions changes in the accretion disk, 
there are two factors, which involves the change delay of the turbulent viscosity. Firstly, 
since the turbulence may be caused by the nonlinear evolution of unstable oscillation 
modes, that for formation of the developed turbulence is required of the characteristic 
time Ti- The first approximation may be thought of as ti is proportional to the build¬ 
up time of instabilities. Secondly, there are delay T 2 , associated with transfer of an 
energy from large-scale pulsations to small-scale pulsations. Hence, as accretion-disk 
parametrs changes (of the temperature and density), value of the actual viscosity time 
delay from instant value of the dynamic viscosity 77 * by characteristic time t — ti + T 2 - 
For standard ct-model AD ~ oiaVlh^, where D is the Keplerian angular velocity. 
The first approximation law of relaxation viscosity 77 to value 77 * may be writen in the 
following form: 

^ ^ rj*-rj . ^ 

dt T ' ^ ’ 

We restrict ourselves to the case of small perturbation with kr ^ 1 and m/r ^ k 
(k is the radial wavenumber and m is azimuthal wavenumber), which allows us to use 





WKB approximation and seek the solution in the form 

/ = fiexp{—iu>t + ikr + irrup) , ( 2 ) 

where u) is the complex frequency of a mode. Equilibrium quantities are denoted by the 
subscript ” 0 ”. 


The influence of the delay of the viscosity. 


In the general case r > 0 we obtain the hve-order dispersion relation. This is 
equation describes hve oscillatory modes. Four of them were considered previously at 
r = 0 [11]. The inclusion of the delay r > 0 gives rise to new oscillatory mode is 
the second-viscous mode, beside Re(c<;) = 0 and Im(u;) < 0 at any values of another 
parameters. Rewriting equation (1) in terms of the result (2) we obtain in the linear 
approximation : 

- = — , (3) 
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For a standard ct-model accretion disk in the radiation-pressure-dominated region 
the increment of all four unstable modes decreases, with the increase of characteristic 
time of the delay r, as indicated by hg. 1. Stabilization of the thermal and viscous 
oscillation modes (Im(c<;) < 0 ) occurs at r > 100 / 0 , while the acoustic mode tend to 
become stable at r ~ I/O. 


The influence of the opacity 

The conditions of development for unstable oscillation modes very strongly depend 
on the model of opacity. We assume that opacity k is a function of a and h (in other 
words, of density and temperature). The linear approximation yelds 
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In the case of the Thomson scattering R — R^s — 0,4sm^/g (Ag- = 0, Ah = 0), 
and for Kramers law R — Rff oc pT~'^ (A^r = 1, Ah = —8). For low temperature 
protoplanetary disks k oc [13] {A^ = 0, A^ = 4). In the ’’hot” limit of Foulkner 

model [14] k cx pT-^/^ {A^ = ’’cold” limit 

k oc (A^ = 1/3, Ah = 59/3). Here /do = Porad/iPorad +Pogas), Pgas IS the gas 

pressure, Prad is the radiation pressure. 

As indicated by hg. 2 the thermal mode of the oscillation become unstable even in 
the case /3o = 0 at A^^ < 0 and A^ > 0. Stabilization of the acoustic-mode occurs at a 
high negative value A^j. However, in the framework of standard model AD [1] acoustic 
modes of the oscillations prove to be unstable both in the inner radiation-dominated 
region (Prad > Pgas, R = Res), and in the external gaseous region {Pgas > Prad, 
R = Rff). 
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The influence of the viscosity 

As the second (elastic) viscosity hq increases, the increment of the sonnd waves 
decreases nntil the imaginary part vanishes at //q = fJ^ocrit- The pertnrbation decay 
(Im(ci;) < 0) when fxo > fiocrit- As shown in the work [11] the qnantity fiocHt only 
weakly depends on /Iq, so that Hocrit = 2 ^ Sz/q- 

The increments and conditions for instability development very strongly depend 
on the parameters 5^ and dh- By now there are many models of tnrbnlent viscosity 
with differents valnes her, dh- For standard o-model Wrip — —ap [1] and = 0, 
5h = 2. For modiheations of o-models = —a{pg/p)^f‘^p {N = const) [2,4,5] and 
, m , 4 + /lo(12-7iV) 

"" 2(l + 3/lo)’ 2(1 + 3/1o) ■ 

The fnnetions Im[ti;(ho-)] and Im[c<;(h/i)] are in qnalitative agreement with Im[c<;(Ao-)] 
and Im[c<;(A^)] accordingly. In the terms of the mentioned models of tnrbnlent viscosity, 
both aconstic modes of the oscillations are nnstable at pq < pocHt, becanse for all these 
models her > 0 and dh > 0. The aconstic oscillations stabilize at large negative valnes 
of da- 

Conclusion 

A linear analysis has been performed to examine the radial-azimnthal instability of 
accretion disks. The main resnlts are as follows: 

First, the existence of nnstable modes is associated completely with the pertnrba¬ 
tion of dynamic viscosity rj = au and, conseqnently, it is determined by the dependence 
of z/(cT, h). If we set fj = 0, all fonr oscillatory modes (inclnding the thermal and viscons 
modes in the radiation-dominated limit) wonld decay with decrement Im((h) r'-j -l/QP. 

Second, the azimnthal pertnrbation wavennmber m in terms of the WKB approx¬ 
imation for a short-wave pertnrbations with kr ^ 1 not affected by the increments of 
nnstable modes and that is responsible for Doppler shift of the freqnency Cj = u — mVt. 

Third, in the radiation-pressnre-dominated accretion disk the thermal and viscons 
un stable modes stabilize at large valnes of the characteristic time scale of the delay r, 
while the aconstic mode tend to become stable at small valnes r. 

Finally, the conditions of development and parameters of the thermal, viscons and 
aconstic instabilities in the geometrically thin and opticaly thick accretion disk are 
strongly depended on the models of viscosity and opacity. In the gas-pressnre-dominated 
accretion disk the thermal mode tend to become un stable with decreases of the valnes 
dfj , Act and with increases of the valnes dh , A^. The aconstic mode at any gas-to- 
radiation pressnre ratio tend to become stable with decreases of the valnes d^j , A^r. An 
additional point to emphasize is that the sonnd waves stabilize at large valnes of the 
second (elastic) viscosity. 
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Figure Captions 

Fig. 1. Dependence of the imaginary part of the freqnency in terms of the angnlar 
velocity Im(c<;)/D on the withont dimensional wavennmber kh^ a),b) and on the withont 
dimensional time of delay rO c),d). 

Fig. 2. Dependence of the imaginary part of the freqnency in terms of the angnlar veloc¬ 
ity Im(a;)/0 of aconstic, thermal and viscons oscillation modes on the valne parameter 
Act a),b) and on the valne parameter c),d) at kho = 1. The solid and short dashed 
lines represent /Sq = 0 and /3o = 1, respectively. 
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